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We calculate production cross sections of a forward quark-gluon pair and of two gluons 
at mid-rapidity in Deep Inelastic Scattering and in high energy proton-nucleus collisions. 
The calculation is performed in the framework of the Color Glass Condensate formalism. 
, We first calculate the cross sections in the quasi-classical approximation, which includes 

■ multiple rescatterings in the target. We then proceed to include the effects of non-linear 

small- x evolution in the production cross sections. It is interesting to note that our result 
for the two-gluon production cross section appears to be in direct violation of AGK cutting 
rules, which is the first example of such violation in QCD. The calculated quark-gluon and 
Oh, gluon-gluon production cross sections can be used to construct theoretical predictions for 

two-particle azimuthal correlations at RHIC and LHC as well as for Deep Inelastic 

Scattering experiments at HERA and eRHIC. 
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I. INTRODUCTION 

Recent prediction of high-p^ suppression in the nuclear modification factor R dA at forward rapidity RHIC 
dAu collisions [1-3] based on the physics of parton saturation/Color Glass Condensate [4-11] has been 
confirmed by the experimental data in [12-16]. The prediction of [2] was based on the calculation of inclusive 
gluon production cross section in DIS and pA collisions. The calculation was first done in the quasi-classical 
framework of McLerran-Venugopalan model including all multiple rescatterings [17] (see also [18-21]). The 
effects of non-linear small- a; evolution [9,10,22] were included in the obtained formula in [23] (see also [24]). 
In [2,3] it was argued that at lower energies/rapidities, where the particle production is given by the quasi- 
classical formula from [17], the nuclear modification factor R pA should exhibit low-pr suppression together 
with a strong enhancement at high-pr, known as Cronin effect [25] (see also [26-29] for similar conclusions). 
However, at higher energies/rapidities, when quantum evolution becomes important, one should expect 
suppression of R pA at all pr [1-3] due to the onset of BFKL anomalous dimension for gluon distributions 
[30]. (It had been earlier suggested in [31] that the forward rapidity region would be most sensitive to small- x 
evolution effects.) Similar argument about enhancement and suppression can be carried through for valence 
quark production cross section calculated in [31]. The suppression has been confirmed experimentally in 
[12,13]. The centrality dependence of the observed suppression was also in agreement with the predictions 
of the Color Glass Condensate formalism [2]. Further developments in the area included an analysis of 
running coupling corrections [32] and a study of similar suppression in di-lepton production [33] (see also 
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[34]). Recently, a more quantitative analyses [35,36] based on the Color Glass Condensate formalism have 
been performed which show good agreement with the data [12,13]. 

Another distinctive prediction of the Color Glass Condensate (CGC) [4-10] is the disappearance of back- 
to-back jets in the low pr < Q s and intermediate pr > Q s transverse momentum regions. While the single 
particle spectra in dAu collisions at RHIC have been successfully described by CGC-inspired models [35,36], 
it is important to go beyond single particle spectra and probe other observables, such as two particle corre- 
lations, in order to map out the region of phase space where CGC is the dominant physics. The inclusive 
two particle (gluon) cross section at high energy is given by the fey-factorization [37] in the high px region 
(pt 3> Q s ) with the gluon distribution function evolving via the BFKL evolution equation [30]. Models 
based on kr factorization have been applied to many different processes, such as the non-flow contribution 
to the elliptic flow observable v 2 in heavy ion collisions [38]. Recently a similar model of two-particle corre- 
lations in dA was used in [39] to predict broadening and disappearance of back-to-back correlations in pA 
(or dA) collisions. The predictions of [39] appear to be confirmed by the preliminary data reported in [40], 
thus strengthening the case for saturation/Color Glass Condensate in dAu data at RHIC. 

Nevertheless, a theoretically rigorous treatment of inclusive two-particle production in DIS and proton- 
nucleus collisions in the low pr region (pt^Qs) has not been performed yet. It is clearly needed in order 
to provide reliable predictions in the px < Q s momentum region, which is the region where the new physics 
of CGC is expected to be most pronounced. (Very recently, there has been a series of articles investigating 
quark-antiquark production in pA collisions using the quasi-classical approximation in the CGC formalism 
[41].) 

Our goal in this work is to derive inclusive two-particle production cross sections using the Color Glass 
Condensate formalism. We start by considering production of two gluons. We assume that the two gluons are 
separated by a large rapidity interval so that their respective rapidities are ordered, yi » y\ . This kinematics 
is, for instance, relevant to the case of two particle production in p(d)A collisions (for example at RHIC or 
LHC) when one of the produced particles is in the mid-rapidity region while the second particle is closer to 
the forward rapidity region. In Sect. II, we derive an expression for two-gluon inclusive cross section in Deep 
Inelastic Scattering (DIS), using the quasi-classical approximation (the McLerran-Venugopalan model) in 
the Color Glass Condensate formalism and making the large- N c approximation to simplify the calculations. 
(The quasi-classical approximation employed here is identical to the one used in [17-19] to describe single 
gluon production.) The result for two-gluon production cross section in the quasi-classical approximation is 
given by Eqs. (1) and (13). We note that a similar expression for two-quark production cross section was 
obtained previously for DIS in [42]. In Sect. Ill, we include the effects of nonlinear small- a: evolution [9] in 
the two-gluon inclusive cross section obtained in Sect. II. The final answer for the two-gluon inclusive cross 
section for DIS is given in Eq. (32). This result can be easily generalized to pA collisions. 

An ansatz for two gluon inclusive cross section including saturation effects was written in [24] inspired 
by fcr-factorization together with AGK cutting rules [43]. We note that the diagrammatic structure of 
our answer in Eq. (32) docs not seem to adhere to AGK cutting rules' expectation for two-gluon inclusive 
cross section. Furthermore, we are unable to cast the expression (32) into fcT-factorized form used in [24]. 
However, the leading twist fey-factorization expression [37,23] can be reproduced exactly from Eq. (32), as 
will be discussed at the end of Sect. III. 

In Sect. IV we calculate the inclusive production of a valence quark and a gluon in pA collisions both in 
the quasi-classical approximation and including the quantum evolution in the target. The rapidities of the 
valence quark and the gluon are assumed to be comparable and large (both quark and gluon are produced in 
the forward rapidity region). The result is given by Eqs. (76- 78). These expressions together with Eq. (32) 
can be used to describe the nuclear modification factor for azimuthal correlations I dAu at any rapidity 
between mid-rapidity and the deuteron beam at RHIC. In particular, Eq. (32) provides the theoretical basis 
for the correlation analysis carried out in [39] . 

II. TWO-GLUON PRODUCTION IN THE QUASI-CLASSICAL APPROXIMATION 

In this Section we are going to derive an expression for inclusive two-gluon production cross section in DIS 
including all multiple rescatterings of the two produced gluons and the quark-antiquark pair on the nuclcons 
in the target nucleus [44,7,8]. A typical diagram contributing to the process is shown in Fig. 1. The two 
produced gluons have transverse momenta k x and fe 2 and rapidities yi and y 2 correspondingly. To simplify 
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the calculations we will consider the case when 2/2^2/1- A more general case of y 2 ~ Ui was considered in 
[45] for two-gluon production at the leading twist level given by fey-factorization. Our goal here is to include 
the saturation effects in the two-gluon production cross section, which means summing all twists. We will 
achieve this difficult task only for a simpler case of 2/2 3> yi, though, in principle, the more general case 
2/2 ~ !/i presents no new conceptual difficulties and is only technically more complicated. 

As shown in Fig. 1 the gluon production process in D1S in the quasi-classical approximation [44,7,8] 
consists of two factorisable stages. First the incoming virtual photon splits into a quark-antiquark pair, 
which emits two gluons in the incoming wave function. (The time scale for this splitting and gluon emissions 
is much longer than the time of interaction with the target.) The whole system multiply rescatters on the 
nuclear target. (In general the gluon emissions can happen after the interaction with the target, as will be 
discussed shortly.) In the quasi-classical approximation considered here the interactions with the nucleons 
are limited to no more than two exchanged gluons per nucleon (see the second reference in [8] ) . Single gluon 
production in the same approximation was calculated for pA collisions in [17] and for DIS in [20]. 

Let us assume that the virtual photon has a large "+" component of the momentum and the nucleus had a 
large "— " component of its momentum. Than the diagram of the process shown in Fig. I is dominant in A + = 
gauge. In this Section we will perform all the calculations in the framework of the light cone perturbation 
theory in A + = gauge [46]. First of all let us explicitly factor out the wave function <3> 7 ^ qq {x Q 1 a) of the 
virtual photon splitting in a quark-antiquark pair of transverse size x Q q = x — Xq with the quark carrying a 
fraction a of the virtual photon's light cone momentum. The wave function <I> 7 ^^(Xqq, a) is a well-known 
function and can be found, for example, in [47,48]. The two-gluon inclusive production cross section can be 
written as 

d a l" A^qqGGX If f 1 * - d<J qqA ^ qqGGX 

¥k~! d Vl cPk 2 dy 2 = 2^ J d2x °" Jo "" ( - o6 ' a) ( Pk 1 dy 1 (Pk 2 dy 2 { ^ ) - (1) 




A. Time-Ordering Rules 



To calculate two-gluon production cross section for a quarkonium scattering on a nucleus, similarly to 
[17,20] one has to consider various possible ordering of the emissions of the two gluons by the qq pair. The 
interaction with the nucleus target can be considered instantaneous compared to long lifetimes of emitted 
gluons. Thus we will denote the moment of interaction with the target by the light cone time r = x + = 0. 
If T\ and t 2 are the times of the emission of the two gluons, the possible emission ordering in the amplitude 
reduces to three cases: (i) both gluons are emitted before the interaction, t\,t 2 < 0; (ii) one gluon is emitted 
before the interaction and the other one is emitted after the interaction, n < < r 2 or t 2 < < t\) (iii) 
both gluons are emitted after the interaction, t\ 1 t 2 > 0. 

The three cases are represented in Fig. 2 for a particular coupling of the two gluons to the qq pair. 
There the dashed line in the middle denotes the (instantaneous) interaction with the target. The dashed line 
comprises a-Uthe multiple rescatterings like the ones shown in Fig. 1. The dotted lines represent intermediate 
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states, which will give energy denominators in light cone perturbation theory [46]. Even though Fig. 2 shows 
a particular way of the gluons' coupling to the qq pair, the conclusions we will draw below about which 
diagrams dominate will be applicable to other couplings of the gluons. 

Let us define the light cone energy of a gluon or a quark line carrying momentum (7c, k + ) as [46] 



E k = k- 



k 2 



2k^ 



(2) 



In Regge kinematics that we consider here the light cone momenta of the gluons are ordered, such that 
k^ 3> k\ and E k2 <C E kl ■ First we consider case (i) in Fig. 2. The diagrams (i)A and (i)B (top and bottom) 
are different only by energy denominators. Therefore, forgetting the rest of the diagram for now, we write 



(i)A 



1 



1 



Ek 2 E kl + E k2 



1 1 

Ek 2 E kl 



(3) 



and 



(i)B 



E kl E kl + E, 



k 2 



(4) 



The intermediate states giving the energy denominators in Eqs. (3) and (4) are shown by dotted lines in 
Fig. 2. Since E k2 <C E kl , Eqs. (3) and (4) imply that (i)A 3> {i)B. Therefore the diagram (i)B can be 
safely neglected. The conclusion we draw from this analysis is that for gluon emissions before the interaction 
(n,T2 < 0) the (longitudinally) harder gluon has to be emitted first, as pictured in the diagram (i)A in 
Fig. 2. 
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FIG. 2. Possible orderings of the emission of the two gluons by the quark-antiquark pair. 



In calculating the diagrams (i)A and (i)B we have neglected the "— " component of the momenta of the 
qq pair, because they are negligibly small. The quark and the antiquark carry a very large "+" component 
of the momentum, of the order of p + 3> k^ ^> k\ , which leads to negligibly small light cone energy E p . 
We have also neglected the change in the "— " component of the target momentum, since the interaction 
with the target took place after the intermediate states which gave the energy denominators in case (i) in 
Fig. 2. This is not the case in the rest of the diagrams in Fig. 2. The "— " momentum/light cone energy of 
the target changes due to the interaction (dashed line). However, since the light cone energy is conserved 
in the final state (t = +oo), the change of the target's "— " momentum is compensated by the change of 
the "— " momentum of the projectile, which is mostly due to appearance of two extra gluons leading to an 
addition of extra E kl + E k2 to the qq wave function's light cone energy. Therefore, the target's light cone 
energy decreases by E kl + E k2 after the interaction. Thus, when calculating the energy denominators of the 
intermediate states after the interaction one has to add the change in the light cone energy of the target to 
the energies of the lines shown in Fig. 2. This is equivalent to subtracting E kl + E k2 in the corresponding 
energy denominators. (This rule is worked out in more detail in Sect. IIIA of [23].) 
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Guided by the rule we just derived we write for the energy denominators of the diagrams in the case (ii) 
in Fig. 2 

{U)A 1 1 = (5) 

E kl E kl - (E kl + Ek 2 ) E kl E k2 ' 

(ii)B 1 — = ——. (6) 

E k2 E k2 — (E kl + E k2 ) E k2 E kl 

As we see from Eqs. (5) and (6) the two diagrams are of the same order, (ii)A <~ (ii)B [<~ (i)A], and 
neither of them can be neglected. (As we will see below diagrams (ii)A and (ii)B are different in the parts 
responsible for the interaction with the target, so while being parametrically of the same order, they are not 
identically equal.) 

Finally, calculating the graphs in the case (iii) of Fig. 2 one arrives at 

*~ ^ E kl + E k2 E kl ElJ ^ 

and 

(iii)B 1 1 - 1 1 (8) 

E kl + E k2 E k2 E kl E k2 

Since E kl 3> E k2 we conclude that (iii) A >C (iii)B. The diagram (iii) A should be neglected. Therefore, we 
derive a rule for late-time emissions, which take place after the interaction (n, ti > 0): the harder gluon has 
to be emitted after the softer gluon. It is interesting to note that this ordering is the exact inverse of the 
ordering giving the leading contribution at early times before the interaction. The rule can also be generalized 
to any number of gluon emissions contributing to the BFKL [30] or, equivalently, dipole evolution [49,9]: in 
the evolution at early times preceding the interaction the gluons are ordered so that the harder gluons are 
emitted before the softer ones [49]. The ordering is reversed for late times following the interaction, where 
the harder gluons should be emitted after the softer gluons to pick up the leading logarithmic contribution. 
This observation was made previously in [23] . 



B. Two-Gluon Inclusive Cross Section in the Quasi-Classical Approximation 

The diagrams contributing to emission of the harder gluon with momentum (k 2 , y 2 ) are shown in Fig. 3. 
(In the following we will refer to this gluon as gluon #2 and to the other (softer) produced gluon as gluon 
#1.) To simplify the color algebra we will continue the calculation in 't Hooft's large- N c limit. Only planar 
diagrams will contribute for gluon emission. Using the notation from Mueller's dipole model [49,50] we denote 
the gluon in the large- N c limit by a double quark line and leave the ends of the gluon line disconnected from 
the quark lines. The latter notation indicates a sum over all possible connections of the gluon to the qq pair. 

As in Fig. 2, the dashed lines in Fig. 3 denote the t = moment of the interaction of the system with the 
target nucleus. However, unlike Fig. 2, in Fig. 3 we depict the squares of the amplitude contributing to the 
total production cross section. Therefore each diagram has two dashed lines corresponding to interaction 
with the target in the amplitude and in the complex conjugate amplitude. The solid vertical lines denote 
the final state at r = +oo. 

Similar to Mueller's dipole model [49] the emitted gluon #2 in Fig. 3A splits the dipole 00 into two color 
dipolcs. The emission of the softer gluon #1 can happen in any of these two dipoles. However, the original 
dipole model [49] was written for the calculation of the total cross sections, where one only has to calculate 
the forward scattering amplitude of the quarkonium. In that quantity all the final state emissions (r > 0) 
cancel, as was shown in [50]. This is not the case for the inclusive production cross section that we want 
to calculate here. All final state emissions has to be taken into account, as shown in Fig. 3. Also, since 
we are interested in gluon production, the momentum of gluon #2 is fixed. Therefore, since we are going 
to perform our calculations in transverse coordinate space, we have to keep the transverse coordinates of 
the gluon #2 different on both sides of the cut. (To obtain the cross section we will afterwards perform a 
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Fourier transform into transverse momentum space.) Thus, the gluon's transverse coordinate is denoted by 
x 2 to the left of the cut and x 2 , to the right of the cut. Then the color "dipole" formed by, say, the lines 2, 
2' and in Fig. 3A would not be literally a dipole since one needs more than two transverse coordinates to 
describe it, but it would still have the color topology of a dipole and we will refer to it as a "dipole" below. 

A B C D 








co 



FIG. 3. All possible emission of the harder gluon #2 by the qq pair. 



Let us start by analyzing the gluon production in Fig. 3A. As was mentioned before, the softer gluon 
#1 can be emitted either off the color "dipole" formed by lines 0, 2 and 2' or off the "dipole" 2, 2', 0. In 
the following analysis we will concentrate on the latter case of emission of gluon #1 in "dipole" 2, 2', 0. (A 
generalization to emission in "dipole" 0, 2 and 2' is straightforward.) We will denote by M (x 2 , x 2 , , x$ ; k x ) 
the cross section of emission of a softer gluon #1 in the "dipole" 2, 2', 0. Than the "dipole" 0, 2, 2' would not 
have gluon emissions in it, but it would still be able to interact with the target. Interactions of the target 
with the line would cancel due to real- virtual cancellations [17,48,23]. Interactions with the lines 2 and 
2' do not cancel: instead they are given by the S-matrix of a 22' quark dipole interacting with the target 
[23,51]. The S-matrix is given by 



S (x 2 ,x 2 ,) = 1 - N (x 2 ,x 2 ,), 
where, in the Glauber-Mueller approximation, the forward scattering amplitude No is [44] 

N (x 2 ,x 2 ,) = l- e - x 22' Q *o ln(V«»/A)/4 



(9) 



(10) 



where x 22 > = \x 2 — x 2 i | and the quark saturation scale in the McLerran-Venugopalan model Q s q [7,8] is given 
by (in the large-A^ limit) 



Q 2 s0 (b) = 2wa 2 sP T(b), 



(11) 

the nuclear profile function 



with p the atomic number density in the nucleus with atomic number A, T( 
with b = (x_ 2 + x 2 ,)/2 and A some infrared cutoff. 

In the diagram in Fig. 3B the softer gluon #1 can not be emitted off gluon 2' in the complex conjugate 
amplitude: that would be suppressed due to the inverse ordering rule we derived in Sect. IIA. Furthermore, 
if the gluon 2' is emitted off the quark line the gluon #1 can not be emitted in the dipole formed by 
lines and 2' due to the same inverse ordering rule. Therefore, if the gluon 2' is emitted off the quark 
line 0, the gluon #1 can only be emitted by the dipole 00 in the complex conjugate amplitude. Therefore, 
the diagram Fig. 3B would bring in a factor of M (x 2 , x , x^; k^ if the gluon #1 is emitted in the lower 
"dipole" . In the same case, in the upper "dipole" , only the dipole 02 would interact with the target bringing 
in a factor of Sq (x , x 2 ) . The diagram in Fig. 3C can be obtained from 3B by horizontal reflection which 
can be accomplished by interchanging x 2 <-> x 2 , . Finally, a similar line of arguments shows that the diagram 
in Fig. 3D brings in a factor of Mo(x ,x ,Xq, k x ) if the gluon is emitted in the lower "dipole". 

Combining the diagrams A-D in Fig. 3 and defining 



TV 



(12) 



we write 



d 2 kx dyi d 2 k 2 dy 2 



(2tt) 3 



d B d x 2 d x 2 > e 





2^-20 




1 2l2'0 


^2'6 ] 




y x 2o 


~4o)' 


{ X 2'0 


_ 4oJ 
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X Mo(Z2,£2',£o;£l)<Sofe2,« 2 ') 



L20 



L20 



20 j 



2 — Mo^XQjXQiki) S (x ,x 2 



"2'0 



P I ' -F M o(£ 0; £2'>£5;Ai) Sb(2o, £2') + =^ • -^■M (x o ,x o ,sco;k 1 ) + (0^6) 



x 2 

x 2'0 



' 20 



20 



X 2 
x 2'0 



(13) 



where B_ = (x + xn)/2 is the impact parameter of the original dipole 00. The term (0 <-> 0) implies that we 
have to add the whole expression again interchanging and to account for the emission of gluon #1 from 
the top "dipole" . 
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FIG. 4. All possible emissions of the softer gluon #1 by the "dipole" 2,2',0. 

Now we have to calculate Mq{x 2 , x 2 > , 2.5; k t ) . To do that let us consider all possible emissions of the gluon 
#1 in the "dipole" 2,2', as shown in Fig. 4. The transverse coordinates of gluon #1 are x_ x and x_ v to the 
left and to the right of the cut correspondingly. 

To calculate all the diagrams in Fig. 4 one has to use the rules of Sect. IIA. Let us illustrate the prescription 
for calculating these graphs for the fairly general case when the gluon #1 is emitted off gluon #2 on both 
sides of the cut. The "dipole" 2,2',0 in Fig. 4A would then split into a "dipole" 1,1', and a quadrupolc 
2, 2', 1, 1'. The interactions of line with the target cancel via real-virtual cancellations [17,48], thus reducing 
the interactions of the "dipole" 1, 1', to the interaction of a real dipole 1, 1' bringing in a factor of So(x 1 ,x 1 ,). 
The interaction of the quadrupole 2, 2', 1, 1' brings in a factor which we will denote Qo{x 2 ,x 2 , ,x x ,Xy). In the 
quasi-classical approximation of McLerran-Venugopalan model this S'-matrix of the quadrupole interaction 
with the target Qo is calculated in the Appendix A yielding (cf. [42]) 

Oo(xa,^,£i,£i0 = e-l^Mi/^iAJ+x^^Mi/^i'AMQjo/d 

| x\ 2 , ln(l/s 2 2'A) + 3-u/ ln(l/xn/A) - x\ v ln(l/s 2 i'A) - x 2 , x ln(l/a; 2 /iA) 
x\ x ln(l/x 2 iA) + x\, v ln(l/x2'i'A) - x\ 2 , ln(l/x 22 /A) - x\ v ln(l/a;ii/A) 



-\x\-y ln(l/x 21 A)+a;^ 1 , la(l/x 2 , 1 ,A)]Q= /4 



-\x\ x , ln(l/x u ,A)+x= 2 , ln(l/x 22 ,A)]Q : 



(14) 



As once can see from Eq. (14), Q (x 2 , x 2 , x x , x_ x ) = 1 which is what one would expect due to real- virtual 
cancellations [17,48,23]. One can also check that 



Qo {2L2 ' —2' i2Lii3Li, 



= p-^22' ln ( 1 /^22'A)Qf /4 



(15) 



corresponding to the S'-matrix of the interaction of the dipole 22' with the target. (Interactions with line 1 
cancel again due to real-virtual cancellations if we put x x = x v .) 

Now let us evaluate the graph in Fig. 4B for the same case of gluon #1 being emitted off gluon #2 on 
both sides of the cut. In the top "dipole" interactions can only take place to the left of the cut giving a 
factor of Sq (x 2 , x x ). In the bottom "dipole" the interactions with line #0 cancel, leaving only the dipole 12' 
to interact with the target which brings in a factor of So (x. 1 ,x 2 /). The diagram in Fig. 4C is evaluated in 
a similar way yielding a factor of So{x 2 i , x lt ) Sq{x_ 2 ,x_ v ). Finally, in the diagram in Fig. 4D only the lower 
dipole can interact with the target giving a factor of Sq(x 2 ,x 2 ,). 
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Combining the factors calculated above for the diagrams A-D in Fig. 4, putting in the contribution of 
gluon emission and summing over all possible connections of gluon #1 to the lines 2,2' and 1 we obtain 



M)(£2,£2',£o,fcl) 



(27T)* 



/ 



d 2 x\ d 2 xy e 



±12 _ ±1'2' 

ry>2 ~2 

x 12 -^l^' 



Qo(x 2 , X_ 2 t , X_i , Xy )S (x 1 ,x v ) 



+So{x 2 ,X 2 ,) - S (£2,£i)<So(£l,£2') ~ Sofe'^lO^ofe^l') 



10 



il'O 



l'O 



Oo (^2 ' ^22' ' —1 ' —1' ) 



X'S'ofe,^!') + S'ofeife') - Qo(£ 2 >£2'>£i>£o) S'ofeijSo) ~ Qo(£ 2 :£2'>£5>£l') -S'ofeo'Sl') 



^l'O 

c 2 - 
''l'O 



Qofe.^'i^i^iO'S'ofe,^!') + S {x 2 ,xq) S (x 2 ,,xq) ~ Qoi^x^,^,^) Soix^XQ) 



-S a {x 2 ,x v ) S a {x 2 ,,x v ) 



ao 



10 



il'2' 



'1'2' 



Q (x 2 , X_ 2 I , X_i , X_y 

) 5 + Sq{x 2 ,Xq) S (x 2 ,,Xq) 



-Soix^Xx) Sq{0C 2 ,,X 1 ) ~ Qo{Xj2,Xj 2l ,Xx s ,X v )S( j {x v ,Xf S ) 

Eqs. (1) and (13), together with Eqs. (16), (14), (9) and (10) give us the two-gluon inclusive production 
cross section for DIS on a nucleus including all the quasi-classical multiple rcscatterings in the large-./V c 
approximation. It is the main result of this Section. 




III. TWO-GLUON PRODUCTION INCLUDING QUANTUM EVOLUTION 



In this Section our goal is to include the effects of non-linear small- a; quantum evolution of [9] into the 
quasi-classical expression (13) for inclusive two-gluon production cross section. We will begin by reviewing 
the non-linear evolution equation and its application to single inclusive gluon production. We will proceed 
by deriving the expression generalizing Eq. (13) by including the non-linear evolution [9] in it. We will 
conclude by verifying that the obtained expression matches onto the standard fey-factorization result [37] at 
the leading twist level. 



A. Brief Review of Small-a- Evolution and Single Gluon Production 



To include the effects of small- a; evolution in the dipole S'-matrix one first defines the S matrix for a quark 
dipolc 00 having rapidity Y with respect to the target as 

S(x j0 ,xq,Y) = 1-N(x ,xq,Y), (17) 

where the forward scattering amplitude has to be determined from the non-linear evolution equation [9] 

N(x ,x 6 ,Y) = N (x ,x,)e-^ H?)* + Ji f dy c -»*. K=?)0--») 

2lT Jo 



I 



d 2 x 2 



oo 



2 2 
X 20 X 20 



W{xja,x 2 ,y) + N(x 2 ,x 6 ,y) - N{x ,x 2 ,y) N(x 2 ,x^,y)}, 



(18) 
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with the initial condition given by Eq. (10) and p being an ultraviolet cutoff [49]. The evolution equation 
(18) resums all powers of leading logarithms of center of mass energy a s Y and all multiple interactions with 
the target, which bring in powers of a 2 A 1 / 3 with A the atomic number of the nucleus [8,17,48,20]. 

Eq. (18) is derived in [9] by resumming a cascade of gluons in the incoming qq wave function, which in 
the large- N c limit turns into a cascade of color dipoles. The emissions are similar to the ones we considered 
in Sect. IIA for the early times preceding the interaction. The difference is that in [9] one resums emissions 
to all orders, without limiting oneself to just two gluons. At the leading logarithmic level the contribution 
of this gluonic (dipole) cascade to the ^-matrix of the dipole-nucleus scattering is given by the solution of 
Eq. (18). However, one can use this cascade to construct other useful observables. 

In [23] it was shown that inclusive gluon production cross section in a dipole-nucleus scattering is given 
by the following formula: 

d a qqA-^qqGX r 

— d 2 k dy — feoo) = J d " B n i( x -o, x -6, Y ': x -i>x 2 ,y)d x\ d x 2 s(x lt x 2 ,k,y), (19) 
where we defined 

( 27T r J ,f~t x ki-£il 2 la -Xj I 2 

x [N G (z 1 ,x j ,y)+N G (^,x i ,y)-N G (z 1 ,z 2 ,y)-N G (x^,x j ,y)] (20) 

in terms of the forward scattering amplitude of the adjoint (gluon) dipole, which in the large iV c limit can 
be easily expressed in terms of the forward amplitude of the fundamental (quark) dipole 

N G {xz,x x ,y) = 2N(x Ol x 11 y)~N 2 {x 0l x 1 ,y). (21) 

In Eq. (19) the quantity ni(x , Xq, Y; x 1; x 2 , y) has the meaning of the probability of finding a dipole 12 
at rapidity y in the original dipole 00 having rapidity Y [49] . It obeys the following equation [49] 

-Y 

«-ifeo;£o> r ;£i'£i>y) = ^feo -x^SfxQ -xi)e~ 2as ln (-^)( y -^) + ^± / dy' e~ 2as H^r)^^ ) 

2vr J y 

r x 2 ~ 

x d 2 X2^^[ni(x ,x 2 ,y';x 1 ,xi,y) + n 1 (x 2 ,XQ,y';x 1 ,x i ,y)}, (22) 

Jp X 20 X 2 q 

which is the linear part of the dipole evolution equation (18) equivalent to the BFKL equation [30]. 

The quantity s(x 1 ,x 3 ,k, y) in Eq. (19) is the cross section for single gluon production by the dipole 12 
scattering on a nucleus at rapidity y with the emitted gluon being the first (hardest) gluon in the gluonic 
(dipole) cascade developed by the incoming dipole 12. Than Eq. (19) has a simple physical meaning: it 
convolutes the probability of finding a dipole in the initial onium wave function which would emit the gluon 
with the probability of the gluon emission by this dipole. 




y y 




FIG. 5. Feynman diagram corresponding to single gluon production cross section given by Eq. (19). Emitted gluon 
is denoted by the cross. 
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To recover the quasi-classical result for single gluon production [17] one has to put Y = y = on the right 
hand side of Eq. (19). That would effectively turn off the quantum evolution giving 



d 2 k dy 



(x o5 ) = J d 2 B s(x ,x 5 ,fc,0) = J d 2 Bd 2 Zl d 2 z 2 e-^<^-^ 



2 

2. i ill, 



Ef_l\i+3 gl -±i . ~ 2 -3 ( e -(2L t -2L ] YQ: a Mi/\x-x ] \K)/2 _ g-^-^fQ^ ln(l/| 2l -^ |A)/2 
I v.. — nr. . |2 I » i" . 1 2 \ 



_ e -(£ 2 -^) 2< 32oln(l/l£ 2 -£jA)/2 + -£ 2 ) 2 Q?o ln(l/| £l -£ 2 |A)/2^ (23) 

which is the quasi-classical gluon production cross section found in [17-20]. 

If the evolution equation (18) is pictured as resumming the so-called "fan" diagrams in dipole-nucleus 
scattering [4,5], than the single inclusive gluon production cross section would correspond to diagrams like 
the one shown in Fig. 5. There the produced gluon, which is denoted by the cross, can be emitted only 
from the top ladder in the diagram. As it turned out, emissions from all other (lower) ladders cancel [23], in 
agreement with expectations of the AGK cutting rules [43] (see also [24]). Thus the evolution between the 
projectile and the produced gluon is just a linear BFKL evolution, as we can see in Eq. (22). The evolution 
between the produced gluon and the target is the full non-linear evolution given by Eq. (18), as can bee seen 
from Eq. (19). 

Before concluding the subsection let us define another useful quantity. Following [49] let 

n2(x ,XQ,Y;x 1 ,x i ,y 1 ,x 2 ,x 2 ,y 2 ) 

be the probability of finding dipolcs 11 and 22 with rapidities y\ and y 2 correspondingly in the original 
dipole 00 having rapidity Y . This quantity obeys the following evolution equation [49] 



( v \ as f Y J -25. in(%&)(Y-y) f .2 

n2{x ,Xo,Y;x 1 ,x 1 ,yi,x 2 ,X2,y2) = -z~ / dye 1 ' A "' a 



X3 2 „2 



'max{yi >!/2 } J P x w x zq 

x [ni(^,x^,y]x 1 ,x i ,yi)ni(x 3 ,o^,y;x 2 ,x.2,y2) +ni(^,^,yix 2 ,x_2,y 2 )ni(x 3 ,^,y-,x 1 ,x 1 ,yi)+ 

+ n 2 (x^,x 3 ,y-,x 1 ,x i ,y 1 ,x 2 ,X2,y 2 ) + n 2 (x 3 ,x 6 ,y;x 1 ,xi,y 1 ,x 2 ,X2,y 2 )] (24) 
which is linear and can be solved after one finds m from Eq. (22). 



B. Two-Gluon Inclusive Cross Section with Quantum Evolution 

Now we have all the essential ingredients necessary to include quantum evolution effects into Eq. (13). 
Similar to the analysis carried out in [23] we will separate all the gluons into the ones which are harder 
(have higher rapidity with respect to the target) than the harder of the two gluons with rapidities y\ and 
y 2 that are going to be produced and into the ones which are softer (have lower rapidity) than the gluon y 2 
(U2 > yi). 

Similarly to the analysis of Sect. Ill A in [23], one can easily conclude that all of the harder gluons can 
be emitted only at early (t < 0) times both in the amplitude and in the complex conjugate amplitude. Due 
to the ordering rule from Sect. IIA of this paper, this implies that these harder gluons have to be emitted 
before gluons #2 and #1. Therefore, we have to distinguish two important cases: 

A. The gluons #2 and #1 are emitted in two different dipoles created by the evolution due to emission 
of gluons which are harder than either gluon #2 or gluon #1. 
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B. The gluon #2 is emitted in a dipole created by the evolution consisting of emissions of harder gluons. 
The gluon #1 is emitted either by one of the "dipoles" adjacent to gluon #2 (as was studied in Fig. 4) 
or in a dipole generated by evolution inside one of these adjacent "dipoles". 

Case A is relatively straightforward. Quantum evolution creating two dipoles of given sizes and rapidities 
at times r < is included in the quantity n% from Eq. (24). Emission of each of the gluons #1 and #2 in two 
independent dipoles is equivalent to the same problem of a single inclusive gluon emission in a dipole-nucleus 
collision as considered in Sect. Ill A and is described by the quantity s from Eq. (20), which also includes all 
the successive evolution generated through emissions of gluons softer than either #1 or #2 [23]. Therefore 
the contribution of case A to double gluon production can be written as 

with B_ = (Xjq + £q)/2 as before. 

Contribution of case B is somewhat more complicated. The probability of finding an early time (r < 0) 
dipole in the original onium in which gluon #2 is emitted is described by the quantity ni from Eq. (22). 
Emission of gluon #2 is then described by the diagrams of Fig. 3 and, equivalently, by Eq. (13). The only 
difference is that now we have to include quantum evolution in the quantity M and in the S- matrix Sq. 
The inclusion of evolution into the S-matrix Sq is accomplished in Eqs. (17) and (18). The inclusion of 
evolution into the quantity Mq requires a separate diagrammatic analysis, shown in Fig. 6. 

Let us first define a quantity M(x 2 ,x 2 , , Xq,Y'i kiTUi) which by analogy with Mq has a physical meaning of 
an inclusive cross section of producing a gluon with transverse momentum ^ and rapidity y\ in the "dipole" 
2, 2', having rapidity Y . To write down an evolution equation for M (x 2 ,x 2 , ,Xq, Y; k^, y\) one has to analyze 
a single step of small- x evolution for this quantity. All the important gluon emissions in the "dipole" 2, 2', 
are shown in Fig. 6. We start the analysis with the diagram in Fig. 6A. Emitting the gluon #4 splits the 
original "dipole" 2,2', into a "dipole" 2, 2', 4 and a dipole 40. Than the gluon #1 can be emitted in the 
"dipole" 2, 2', 4, which would bring in a factor of 

M(x 3 ,gs 3 ,,:c 4 ,y;k 1 ,yi) 

with y the rapidity of gluon #4. In this case all interactions in dipole 40 cancel. Alternatively the gluon #1 
can be emitted in the dipole 40, which would bring in the familiar factor of 

J d 2 x a d 2 x b nx(^^y\^,^,y\)s(^,^k 1 ,yx) 

from Eq. (19) describing a single inclusive gluon production in a dipole-nucleus scattering. In this second 
case interactions of the "dipole" 2, 2', 4 with the target would not completely cancel. Instead they would 
bring in a factor of S(x 2 ,x 2 ,, y) = 1 — N(x 2 , x 2 ,, y) corresponding to interaction of dipole 22' with the target. 



ABC 




FIG. 6. Diagrams describing one step of evolution for M. 
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In the diagram shown in Fig. 6B the interaction of gluon #4 with the line 6 to the right of the cut gets 
canceled by the diagram similar to the one in Fig. 6E but with the gluon #4 connecting to the line instead 
of line 2 to the right of the cut [50]. Therefore, gluon #4 interacts only with line 2' on the right hand side of 
Fig. 6B and only with line 2 on the right hand side of Fig. 6E. Due to the inverse ordering rule for late time 
emissions from Sect. IIA the gluon #1 can not be emitted in the dipole 24 in graph B of Fig. 6. (In Fig. 6E 
the dipole 24 is not present in the final state and can not contribute to gluon production.) Therefore, in 
both graphs B and E the dipole 24 brings in only a factor of S(x 2 ,x 4 , y) = 1 — N(x 2 ,x 4 , y) into the evolution 
equation we are constructing. The other "dipole" 4, 2',0 would than bring in a factor of 

M(x 4 , £2 ,, £6 ,y; fc^). 

Diagrams C and F can be obtained from B and E by interchanging 2 <-► 2'. Finally, in the diagram D 
interactions of gluon #4 with the line cancel due to real-virtual cancellations of [50] that we have just 
employed in graphs B and E. The remaining interaction with lines 2 and 2' shown in Fig. 6D does not split 
the "dipole" 2, 2', 0. Due to inverse ordering rule of Sect. IIA there will be no softer gluon emissions in 
"dipole" 2, 2'4, such that all subsequent evolution will take place only in "dipole" 2, 2', 0. Therefore, diagram 
D contributes only to virtual corrections, along with the usual virtual corrections at r < in dipoles 20 and 
2'0. 

Combining the contributions of all diagrams in Fig. 6 we write the following evolution equation: 



M(x 2 ,x 2 , ,£o> y ;^i>2/i) = e 



v ^ ' d(x 2 ,x 2 ,,X£,k 1 ,y- l ) + 



~2ir 



d X4 dye 



~a s In 



2 - ) (Y-y) 




M( £2 , S2 ,, & „; il ,,0 + /AA»,fe, S ,^, a ,» I )»fc, 3 ,i,, W )(l-«fe, & ,ri] 



^42 

C 42 



40 , 



.( 



±42' 



\ r 2 

\ x 42' 



' 12 



Af(£4,^2'.^o,y;^i,yi) [i - ^(x^^y)] - 



^42 _ £*A. 245 \ M (x 2 , £4 , a , y; kl , yi ) [1 - AT(z 2 ,,2U,2/)]). 

X 42 X 42 , J y#42' X 40 / ) 



(26) 



The only thing left to do to complete our analysis is to determine the initial condition for the evolution 
equation (26), which we denoted d(x 2 ,x 2 ,,x -,k 1 ,yi). This quantity is the gluon production cross section in 
the scattering of a "dipole" 2, 2', on a nucleus with the small- x quantum evolution included, in which the 
emitted gluon (ki,yi) is the first (hardest) gluon in the gluonic cascade resumming the quantum evolution 
of Eq. (18). Since the emission diagrams are the same as in Fig. 4, the quantity d(x 2 ,x 2 ,,x Q -,k 1 , y\) should 
be given by the expression similar to Eq. (16), where the dipole and quadrupole S'-matrices So and Q have 
to be replaced by their evolved values. We, therefore, write 



d{±\2, , 2±®i k.ii yi) 



(2tt) s 



d z x x d z xy e *£r£n' 
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2Li'2' 
x 2 



Q(lL2i 3L2't3Lii ±\vi Vi) S{x_-y,x_y, yi) 



+S(x 2 ,x 2 ,,yi) - S(x 2 , x x , yi ) S{xi , x 2 , , y 1 ) - S(x 2 ,, x v , yi) S(x 2 , x v , yi) 



2Lw ^ 
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l'O 



l'O 



Q (2L2 ! 2L2' ) S-i ) 3Lv 1 2/1 ) 



xS(xi,x v ,yi) + S(x 2 ,x 2 ,,yi) - Q(x 2 , x 2 ,, x 1; x$, yi) 5(x 1; xq, yi) - Q(x 2 , x 2 ,, xq,x v , yi)S(xjQ,x v , yi) 
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il2 ifc.1'0 



l'O 



Q(^2 1^2') ^li^i'i S(x_i, x v , yi) + S(x 2 ,Xq, yi) S(x 2 , ,XQ,yi) — Q(x 2 ,x 2 , ,Xi, Xq, yi) 



x %n?6iJ/i) - S(x 2 ,x v ,yi) S(x 2 ,,x v ,yi) 
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~2 

X 10 X 1'2' 



Q(Sl2j £2'> -Ei) •£].' j 2/1 ) ^(-Eij -Ei' ! ?/i) 



+5 , fe>^6-yi) 5 'fe'.^0'yi) - S^x^y-i) S(x 2 , ^x^yx) - Q(x 2 , x 2 ,, x$, x v , yi) S{x v , xq, yi) 

Indeed S(x 1 ,x 1 ,,yi) in Eq. (27) is given by Eqs. (17) and (18). The reason why inclusion of evolution 
just corresponds to replacing the Glauber-Mueller expression (9) for S by the fully evolved Eq. (17) has 
been discussed before in [23]. It was observed there that real- virtual cancellations for the gluon emissions 
contributing to the dipole evolution discussed in [50] act very much like the real-virtual cancellations for 
Glauber-Mueller multiple rescatterings [17]. Namely, if interactions of exchanged Coulomb gluons with a 
given quark line cancel in the multiple rescattering (Glauber-Mueller) picture, than emissions of an s-channel 
gluon by the same quark line at early and late times on both sides of the cut would also cancel [50] . One 
can show that interactions with quark and gluon lines that contribute in multiple rescattering case would 
also contribute in the case of evolution. In the end one concludes that inclusion of quantum evolution can 
be accomplished by replacing So from Eqs. (9) and (10) by S from Eqs. (17) and (18) [23]. 

To calculate Q(x 2l x 2 i,x^, Xy, y\) we have to write down an evolution equation for the S-matrix of the 
evolution of quadrupolc 2, 2', 1, 1'. An evolution equation involving a color quadrupole was derived before 
in [50] to reproduce the BJKP equation [52] for four reggeons in the framework of the dipole model [49]. 
The equation derived in [50] corresponds to off-forward evolution for dipoles in the presence of a single 
quadrupole, with all the evolution included in the dipoles. Therefore it should not be compared to the 
equation we are about to write down, since we are interested in the evolution inside the quadrupole. 

We are going to derive an evolution equation for the quadrupole S-matrix Q(x 2 , x 2 , , x l7 x v , yi) including 
all the non- linear evolution effects. The initial condition for evolution of Qix_ 2 ,x_ 2 , ,x_ x ,x_ v ,y\) is given by 
Eq. (14). To write down one step of the evolution we first redraw the quadrupole as shown in Fig. 7. Instead 
of the amplitude squared pictured on the left of Fig. 7 we will use a form of quadrupole cross section similar 
to the forward amplitude shown on the right of Fig. 7. Obviously the diagram on the right of Fig. 7 preserves 
the color structure of the quadrupole. All the interactions with the target in the graph of the left happen 
along the dashed lines at time t — 0. On the right of Fig. 7 we merge two dashed lines from the graph on the 
left into one dashed line with interactions in it. This way, a real interaction with a nucleon where a single 
gluon is exchanged in each of the dashed lines in the left graph of Fig. 7 becomes a two-gluon exchange 
(diffractive) interaction in the dashed line in the graph on the right of Fig. 7. Again the picture is similar 
to the forward amplitude calculation. 




1=0 

FIG. 7. Redrawing the quadrupole interaction amplitude in the form convenient for including quantum evolution 
(see text). 

One step of the quadrupole evolution in the representation of Fig. 7 is shown in Fig. 8. The step consists 
of an emission of a single gluon #3. For instance, in Fig. 8A the gluon #3 splits the original quadrupole 
2, 2', 1, 1' into a quadrupole 3, 2', 1, 1' and a dipole 23. Fig. 8B gives a similar contribution. In Figs. 8C-F we 
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drew the gluon #3 line as disconnected at dashed line: the gluon line is indeed implied to be connected and 
continuous. The disconnected line is just a notation which we borrowed from [50] for, say, a gluon emitted 
in dipole 21 and absorbed in dipolc 11' in Fig. 8C. However, when connecting the two parts of the gluon 
#3 line in Figs. 8C-F one has to be careful to pick up the leading \arge-N c contribution. For instance, 
in Fig. 8C the leading term consists of gluon #3 splitting the original quadrupole 2,2', 1, 1' into a dipole 
31 and a quadrupole 2, 3, 1', 2'. The dominant contribution of Fig. 8E consists of gluon #3 splitting the 
quadrupole 2, 2', 1, 1' into "dipoles" 2, 2', 3 and 1, 1', 3, in which interactions with line 3 cancel due to real- 
virtual cancellations. Leading terms for other graphs in Fig. 8 can be obtained in similar ways. Combining 
them with virtual corrections yields 

Q( x 2, x 2', x li x l',yi) = e' as "V ? )V1 Qo(x 2 ,X 2 ,,X 1 ,X 1 ,) + 
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(28) 

Note that due to real-virtual cancellations the last two terms in Eq. (28) corresponding to diagrams E and 
F in Fig. 8 contain only the dipole S'-matrices in them. The relative signs of various terms in Eq. (28) are 
easier to determine in the representation of the interaction on the left hand side of Fig. 7 keeping in mind 
that gluon emissions at r < and r > come in with different signs. 









FIG. 8. One step of the quadrupole evolution. 
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As we have already mentioned the linearized version of Eq. (28) (see Eq. (38) below) should be combined 
with Eq. (49) in [50] to complete the description of the BJKP evolution in the framework of the dipolc 
model. (By linearizing Eq. (28) we mean substituting S = 1 — N in it (see Eq. (17)) while keeping only 
terms linear in Q and in N, i.e., neglecting products like QN and NN.) However, as was discussed in 
[50], the contribution of linearized Eq. (28) to BJKP evolution is likely to be small: as we will see in the 
next Subsection, the linearized version of Eq. (28) is almost equivalent to the BFKL equation. Therefore, 
its solution is likely to grow with energy just like a single BFKL pomeron, which is much slower than the 
double BFKL pomeron exchange. On the other hand, a solution of the full BJKP evolution equation for 
four reggeons presented in [53] has an intercept much smaller than that of a single BFKL pomeron exchange 
making the contribution of Eq. (28) potentially important for BJKP evolution in the dipole model. 

Let us check Eq. (28) for consistency with our earlier results. First note that if x 2 = x 2 , the interactions 
with line 2/2' would cancel and the following equality should be true 



1 ' —1' 5 



2/1 )■ 



(29) 



^From Eqs. (15) and (10) we can see that Eq. (29) is certainly true for the initial conditions for Eq. (28) 
given by Qo from Eq. (14). Now, as one can explicitly check, putting x 2 = x 2 , in Eq. (28) (with x 22 > — > p in 
the exponent) and assuming that Eq. (29) is true one readily recovers Eq. (18). Thus Eq. (28) consistently 
maps onto Eq. (18) in the limit of Eq. (29). 

Using Eq. (29) in Eq. (27) and remembering that due to Eq. (21) the S'-matrix of a gluon dipole Sq can 
be expressed in terms of the S'-matrix of the quark dipole as 



Sg{x ,x 1iV ) = S' 2 (x ,x 1 ,y) 



(30) 



we observe that 



d(xjQ } x q, x i, k , y) — s{Xf ) , x_i, k,y), 



(31) 



which verifies that Eq. (27) is consistent with Eq. (20). 

Eq. (28), when solved to find Q, can be used to construct d in Eq. (27), which can then be used as initial 
condition to the evolution equation (26). The quantity M(x 2 , x 2 i , xq,Y; k x , y\) obtained this way can be 
used in Eq. (13) instead of Mq along with S from Eq. (17) to obtain the evolved contribution of case B 
considered above to the two-gluon production cross section. Together with the contribution of the case A 
from Eq. (25) it gives us the following expression for the double gluon production cross section in a quark 
dipole-nucleus scattering [with B_= (x + x$)/2] 
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(32) 



Eq. (32) is the central result of this Section. Together with Eqs. (1), (28), (27), (26), (22), (24) and (18) 
it gives us the expression for the two-gluon inclusive cross section for DIS on a nucleus with the effects of 
nonlinear evolution (18) included. 
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FIG. 9. Feynman diagrams corresponding to double gluon production cross section given by Eq. (32). Emitted 
gluons are denoted by crosses. 

The structure of Eq. (32) is illustrated in Fig. 9, where, if one pictures the evolution of Eq. (18) as 
resumming fan diagrams, the diagrams correspond to the first (case A) and the second (cases B and C) 
terms in Eq. (32). The first term in Eq. (32) corresponds to splitting of the original linear evolution in two, 
which is described by Eq. (24) for n 2 - Then each of the two ladders independently produces a gluon with 
all the possible splittings happening afterwards. This is illustrated in Fig. 9A. The second term in Eq. (32) 
corresponds to nonlinear evolution successively producing both gluons, after which all possible splittings are 
allowed, as shown in Figs. 9B and 9C, where we have divided the nonlinear evolution into the linear (Fig. 
9B) and nonlinear (Fig. 9C) parts. The linear evolution in this second term in Eq. (32) is given by m from 
Eq. (22) and by M from the linear part of Eq. (26). This linear evolution leads to production of both gluons 
#2 and #1 and is illustrated in Fig. 9B. The initial conditions for Eq. (26) are nonlinear, given by Eqs. (27) 
and (28). They include ladder splittings and are pictured by the fan diagram in the lower part of Fig. 9B. 
One should note, however, that Eq. (26), while being linear in M, has extra factors of 1 — TV on its right 
hand side. That means that evolution of M includes ladder splittings between the gluons #2 and #1, one 
of which is shown in Fig. 9C. There the evolution leading to creation of gluon #2 is still linear since it is 
still given by n\ in the second term on the right hand side of Eq. (32). However, since the evolution in the 
rapidity interval between the emitted gluons (evolution of M) is nonlinear, splittings are allowed between 
the gluons #2 and #1, as depicted in Fig. 9C. 

The diagrams A and B in Fig. 9 are the same as would have been expected from AGK cutting rules [43] (see 
also [23] and [54] for similar correspondence between the dipolc model results and AGK rules expectations). 
However, the diagram C in Fig. 9, while being included in Eq. (32), is prohibited by AGK cutting rules. 
Therefore, we seem to observe a direct violation of the AGK rules in QCD. Since AGK rules have never been 
proven for QCD, one should not be too surprised that they do not work here. It is interesting to note that 
AGK violation sets in at the level of the 2-gluon production: single gluon inclusive production cross section 
calculated in [23] adheres to AGK rules and so does the diffractive DIS cross section calculated in [54] . 

The violation of AGK cutting rules in Eq. (32) is due to non- linear terms in Eq. (26), which are in turn due 
to late time (after the interaction) gluon emissions at light cone times t > 0. These terms were not important 
for the calculation of the total cross section in the dipole model [49] : there they were found to cancel [50] . 
Thus if one would try to construct an analogy between the fan diagrams [4] and dipole calculations [9] based 
on the correspondence of total cross sections, one would omit such terms. Since the fan diagrams seem to 
adhere to AGK rules, this omission would lead to the erroneous conclusion that AGK rules should work for 
the production cross section. However, as we have seen above, these late time emissions are important for 
single [23] and double inclusive gluon production, violating the AGK rules for the latter. What appears to 
fail here is the one-to-one correspondence between the fan diagrams and dipole calculations. 

Another difference between our result (32) and the direct application of AGK rules to calculating inclusive 
cross section done in [24] is that nonlinear splittings may start exactly at the point in rapidity when the 
softer of the produced gluons is emitted in the diagram B or exactly at the point of emission of both gluons 
#1 and #2 in the diagram A. Similar discrepancy was already observed when comparing the single gluon 
inclusive cross section calculated in [23] to the results of [24]. 

In comparing our result with the formula obtained in [24] (see Eqs. (30) and (36) there) we note that in a 
general case we could not cast Eq. (32) in the A;T-factorized form of [24]. Again, this distinguishes the case 
of two-gluon production considered here from the case of single gluon production from [23] . 

Indeed our result, given by Eq. (32), is rather complicated, especially keeping in mind that one has to 
first solve evolution equations (22), (24), (18), (28), and (26) in order to obtain the desired two-particle 
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production cross section. In order to make Eq. (32) easier to implement, it is highly desirable to find some 
way of simplifying it. Unfortunately, we could not find any simplification of Eq. (32) in the general case. 
Nevertheless, in certain kinematic regimes Eq. (32) may be simplified. For instance, if the center of mass 
energy of the collision is not too high or if the transverse momenta of the produced gluons are sufficiently 
large (|fci|, \k 2 \>Q s ), the nonlinear saturation effects, such as ladder splittings, could be neglected. This 
implies that the diagrams in Figs. 9 A and C are small with the linear part of the diagram B dominating 
the cross section. This is the well-known leading twist result [37], which we will derive from our Eq. (32) in 
Sect. IIIC below. 

In the opposite kinematic regime of very large center of mass energy of the collision and not too high 
gluons' transverse momenta saturation effects become important. There one can note that in Figs. 9 B and 
C the evolution between the projectile and the (harder) gluon #2 is linear and is given by a single BFKL 
ladder exchange. On the other hand one can show [24] that the diagram in Fig. 9A is dominated by the 
contribution where the triple pomeron vertex in the evolution between the projectile and the gluon #2 is 
all the way up at the projectile's rapidity. Therefore, the evolution between the projectile and the gluon 
#2 in Fig. 9A is given by a double pomeron exchange, and is thus energetically more favorable than the 
single pomeron exchange of Figs. 9 B and C. Since the rest of the three diagrams is parametrically the 
same, one concludes that Fig. 9A dominates in this regime, as was originally shown in [24]. Keeping only the 
corresponding first term on the right hand side of Eq. (32) would significantly simplify the calculation of the 
cross section: since an analytical solution of Eq. (24) for n 2 exists [49] , one would only need to find a solution 
of Eq. (18), for which there is a number of analytical and numerical results in the literature. However, one 
has to be careful in neglecting the diagrams in Figs. 9 B and C. If one is interested in azimuthal two-particle 
correlations of the produced gluons, than the contributions of graphs in Figs. 9 B and C might be more 
important than the contribution of Fig. 9A even deep inside the saturation region [39]. 



C. Recovering the Leading Twist Result 

Let us show that Eq. (32) reduces to the usual "leading-twist" fc^-factorization result [37] in the limit of 
large transverse momenta of the produced gluons. Large transverse momenta correspond to small transverse 
distances. For short transverse distances all the evolution equations written above should be linearized since 
all the non-linearities would be negligibly small. Therefore we can right away neglect the first term on the 
right hand side of Eq. (32), which contains a splitting (which is a non-linearity) in the evolution between 
the target and emitted gluons as shown in Fig. 9A. In the remaining second term on the right of Eq. (32) 
we can put all S « 1 to obtain the linearized expression 
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(33) 



To linearize the evolution equation for M (Eq. (26)) we start by linearizing its initial condition given by 
d. It is determined by Eq. (27) with Q given by Eq. (28). Therefore, we start with the initial conditions 
for Eq. (28) given by Eq. (14). Expanding Eq. (14) to the lowest order in the transverse separations (or, 
equivalently, in Q s0 ) we obtain 
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Qo(x 2 , x_2' i Hi ? x i> 



1 - n (x 21 ) - n (x 2 , v ) - n (x 22 ,) - n (x lv ) + n {x 2V ) + n (a; 2a ), (34) 



LO 



where the two-gluon exchange amplitude n is given by the first term in the expansion of N from Eq. (10) 

(35) 



M x 2i) = \ x liQlo ln ~\"- 



Let us assume that Eq. (28) independently includes linear BFKL evolution in each of the n 's in Eq. (34), 
such that the fully evolved linearized quadrupole amplitude Q is given by 
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with n determined by the linearized version of Eq. (18) corresponding to BFKL evolution [30] 
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After lengthy algebra one can show that q(x 2 , x 2 , , x x , x v ,Y) from Eq. (36) satisfies the linearized version of 
Eq. (28) 
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(38) 



This proves that the ansatz of Eq. (36) is indeed the correct linearized quadrupole amplitude Q. Therefore, to 
construct the initial conditions for the linearized version of Eq. (26) we should substitute q{x_ 2 , x_ 2 , , x_ x , x v , Y) 
from Eq. (36) into the linearized version of Eq. (27) obtaining 
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-n(x 2 ,x 2 ,,yi) - nix^Xy^x) 



+ 



^10 _ ±1'0 
2 * 2 

x w x i'6 



±12 



'12 



^l'O 
X 1'0 



"(^2,^i',2/i) sg.j/i) - n(x 2 ,x 5 ,yi) - nfe,^,,^) 



±10 ±1'2' 



2 2 
• T 10 X V2' 



n(^i,s 2 '>yi) +«(^i',«o.j/i) -"fe';«0'2/i) - 



}• 



(39) 



Here we consider scattering on a large nucleus. The amplitude nix^Xy^y) depends on the transverse size 
of the quark-antiquark pair xiv as well as on the transverse position of the dipolc. However, the transverse 
position of the dipolc is given by the overall impact parameter B_ in the scattering process. The integrals over 
x_ x and x v in Eq. (39), while formally going out to infinity in the transverse direction, are indeed effectively 
limited to the typical hadronic size on which the concept of a gluon still makes sense. The dependence of n 
on B_ is smooth for a large nucleus, slowly varying on transverse distances of the order of the typical hadronic 
size. Therefore we write (see [23] for a similar approximation) 



n{x x ,x 2 ,y) « n(x 12 ,B,y). 
With the help of Eq. (40) we rewrite Eq. (39) as 
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(41) 



where V 2 , is the transverse coordinate gradient squared. Eq. (41) is the initial condition for the linearized 
version of Eq. (26). The latter can be obtained from Eq. (26) by putting all N = in it, which yields 
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where we have introduced a linearized amplitude M denoted by 
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The form of Eq. (41) provides us with the following ansatz for the solution of Eq. (42): 

m{x 2 ,x 2 , ^^-^^yi) = f{x 3 a,Y;k 1 ,y 1 ) + /(x 2 , 5 , Y; fc 1; yi) - f(x 22 ,, Y; k lt y x ) 



(42) 
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with / some unknown functions. Substituting the ansatz of Eq. (44) into Eq. (42) one can see that it is a 
solution of Eq. (42) iff 

/(£ 2 i>^;£i.yi) = \ J d 2 x a d 2 x b n 1 (x 2 ,x 1 ,Y;^,x b ,yi)s(gc a ,x b ,k 1 ,yi). (45) 
The final answer for m is 

™fe 2 :£2'>S)> y ;^i'yi) = \ J d 2 x a d 2 x b [ni(x 2 ,x 6 ,F;x Q ,£ fc ,yi) + ni(x 2 ,,x J o,Y;x a ,x b ,yi) 

-ni(£ 2 ,£ 2 ,,y;£ a ,£ 6 , yi)} s(^ a ,^ 6 , k^yi). (46) 

Similar to Eq. (40) we rewrite 

"lfe,^',*'; 2/i) — > m(X22'>^,^;^o,^6,2/l) (47) 

with the impact parameter B_. Again we assume that for a large nucleus n\ is a slowly varying function of 
B. Using Eq. (46) as M in Eq. (33) than yields 
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feoo) 
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X _L I ^ e -i«b-£in^Ji^MJi±M^ ^AV^ 1 (z,B, K ;i a ,i i , !fl ) S (i a ,i t ,* 1 , J1 ). (48) 
Similarly we rewrite Eq. (20) as 

s{x a ^k^yi) = (|^7J / d 2 «,e-^ln( |ffl ~^ + £atl ) V 2 W N G (w, B, Vl ). (49) 
In the linear regime 

N G (w,B, yi ) « 2nfeB,)/ 1 ) (50) 
with n taken from Eq. (37). Therefore, linearized s can be obtained from Eq. (49) using Eq. (50) 



siin^Sk^Vi) = ^^2 / d 2 we-^ In ""f^f + ~ J ) V 2 w n(w,B, yi ). (51) 



Defining b ab = (x a + x h )/2 we relabel the variables of m in Eq. (48) as [49,23] 

m(z,B,y 2 ;x a ,x b ,y 1 ) — > n 1 (z,x ab ,B-b ab ,y 2 -y 1 ). (52) 
Than the integrals over x a and x b in Eq. (48) can be written as 

J d 2 x ab d 2 b ab V 2 Z ni (z, x ab ,B_ — b ab ,y 2 -yi) su„ {x a ,x b ,k 1 ,y 1 ) 

= J d 2 x ab d 2 b ab \7ln 1 (z,x^ b ,b ab ,y 2 -yi)si m (x^,x b ,k 1 ,y 1 ), (53) 

where we have put the index indicating that we have to use a linearized amplitude s from Eq. (51) and 
shifted b ab by B. Using the explicit solution of Eq. (22) [49] 
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with the eigenvalue of the BFKL equation [30,49] 

X (A) = V(l) - ^(A) - 1^(1 - A) (55) 
and with the help of Eq. (51) we rewrite Eq. (53) as 

^2 V z J d2w ^ d2h ab^n 1 {z,w,b ab ,y 2 - yi)^j L kl (w)n(w,B, yi ), (56) 

where we have defined the operator for Lipatov's effective vertex [30,24] 

L k (z) e^V! e~ ik -- z - v! • (57) 



Performing the integrations over x x and x v in Eq. (48) in a similar manner we finally obtain 

^ y d 2 Bd 2 zd 2 w d 2 & 2 ^ni(x o5 ,z,6 2 ,r-y2)J £fc 2 fe) 
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d 2 6i ^nifz, w,b 1 ,y 2 - yi) ) L kl (w)n(w,B, yi ). (58) 

Eq. (58) has the structure of three BFKL ladders (two factors of n\ and one factor of n) with two Lipatov 
vertices, which are responsible for production of gluons, inserted between them. It is equivalent to the 
fcr-factorization prediction for two-gluon production from a single BFKL ladder [see Eq. (30) in [24]]. We 
have, therefore, proven that at the leading twist (large kr) level, our two-gluon inclusive production cross 
section (32) reduces to the conventional fcr-factorized expression (58). 

Two comments are in order here. First of all, it is a little worrisome that in order to recover the conventional 
leading twist result of Eq. (58) we had to expand the initial conditions for the evolution of quadrupolc 
amplitude Q given by Eq. (14) to the lowest order, as shown in Eq. (34). Of course, by doing so, we have 
shown that the leading twist formula (58) is included in our full expression (32). Nevertheless, taking a 
solution of an evolution equation at very short transverse distances, does not necessarily imply doing the 
same to initial conditions of evolution. For instance, if we are interested in the dipole amplitude N(x J0 , x_q, Y) 
at small Xqq we have to solve the linear part of Eq. (18) with the full initial conditions given by Eq. (10) and 
not with the leading order initial conditions given by Eq. (35). The kernel of Eq. (18) involves integration over 
all transverse sizes including large sizes where multiple rescatterings are important and have to be included. 
Multiple rescatterings become important at lower energies than the small- a: evolution and thus have to 
be included as initial condition even for linear (BFKL) evolution equation. (As one can show, multiple 
rescatterings in the quasi-classical limit become important at rapidity y m uit ~ lnl/a s , while the BFKL 
evolution becomes important at ysFKL ~ l/ a s-) Thus, expanding the initial conditions of Eq. (10) would 
not be justified if one is interested in the small x g of the amplitude N(x JQ , x$, Y). The effects of saturation 
in the initial conditions on the short distance/large fcx behavior of the amplitudes and gluon production 
cross sections have been studied before in [2,1]. It is exactly these effects that bring in suppression of the 
nuclear modification factor R pA for the gluon production [1-3]. Therefore, taking the large ki and k 2 limit of 
Eq. (32) more carefully may result in an expression different from the fey-factorization formula of Eq. (58). 

The second observation one has to make is that by now the reader can appreciate the tremendous simpli- 
fications one needs to make [e.g. Eq. (34), linearization of all evolution equations, etc.] in order to recover 
the fey-factorization formula of Eq. (58). This is strikingly different from the case of single inclusive gluon 
production considered in [23] . There the obtained expression for the cross section was cast in a fcy-factorized 
form without making any linearization assumptions, i.e., without taking the leading twist (high-fcx) limit. 
This /cr-factorization result of [23] was, indeed, unexpected and very puzzling. However it is also interesting 
to observe that it does not hold for the double inclusive gluon production cross section (32). This leaves 
us guessing whether the preservation of fc^-factorization in the formula from [23] for the single gluon in- 
clusive production cross section after multiple rescatterings and small-x evolution had been included (see 
Eq. (19) above) was just incidental. A similar breakdown of fey-factorization has been observed recently for 
qq production in pA collisions [41]. 
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IV. VALENCE QUARK-GLUON PRODUCTION IN PROTON-NUCLEUS COLLISIONS 



In this section we calculate the cross section for production of a valence quark and a gluon in high energy 
proton nucleus collisions. Both produced quark and gluon assumed to have similar rapidity and to be in 
the proton (deuteron) fragmentation region. In this case, one can treat the proton (deuteron) as a dilute 
system of partons while the target nucleus is treated as a Color Glass Condensate. The produced quark and 
gluon then fragment into jets which can be measured. Previously, this approach has been used to calculate 
valence quark, photon and di-lepton cross sections in proton (deuteron)-nucleus collisions [31,33]. Here, we 
extend this formalism to production both of a quark and a gluon. Explicitly, we calculate the differential 
cross section for the following process: 



q(p) A ^ q(q) g{k) X 

given by the amplitude 

M(q,k;p) =< q(q) g(k) out \q(p) in >=< Qout\a ut{k)b out (q)b\ n {p)\Q m > 



(59) 



(60) 



which, using the LSZ reduction formalism, can be written as (we set the renormalization factors equal to 1 
since we are working at the leading order in a s ) 

M=g [ d 4 x d 4 y d 4 z d 4 r d 4 f e^ q - z+k - r - p - y) u(q)[i z ]S F (z, x^t S F (x,y)[i y }u(p) 



Gl» p (x,f)D^(f,r)e,(k) (61) 

where Sf, G vp are the quark and gluon propagators in the classical field background and D PIL is defined 
such that 



J d 4 rGl^{x,r)DZ{r,y)^8l5^5\x-y) 



(62) 



where G° p is the free gluon propagator. This amplitude is shown in Fig. 10 where the quark and gluon lines 
with a thick dot represent the propagators in the background field as illustrated in Fig. 11. 




FIG. 10. Production of a quark and a gluon including multiple scattering from the target. 



I 




FIG. 11. Multiple scattering of a quark or gluon on a target. 
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To proceed further, we write the the propagators in the above amplitude in momentum space. The 
amplitude is 



M=9 J (t^ 4 (to 4 ° (g) * ^fa fc 07^ c ^(fc 2 ,p) ? u{p)G%{k 2 - fc 1 ,fc 3 )^(fc 3 ,fc)e M (fc) (63) 

The quark and gluon propagators in the classical background field are already known [55,56]. It is useful 
to separate the free and interacting parts of the propagator in the following. Therefore, we Define the 
interaction part of the propagators in momentum space as (and suppressing the color factor for the moment) 

S F (q,p) EE (27r) 4 5 4 (p - q) S° F (P) + S° F (q) Tf(q,p) S° F (p) 
G^(q,p) ee (27r) 4 5 4 (p - q) G°^(p) + G°f(q) r g (q,p) G°" v (p) (64) 

where the free propagators are 



<4f and G %( k ) = T 
p z p k 



. f)iik v + ^vkn 
9tw + 



rj ■ k 



(65) 



and r]^ is the light cone gauge vector so that rj ■ A = A~ = defines the gauge we are working in. In 
this gauge the interaction part of the gluon propagator in Eq. (64), denoted here by T g (q,p), is diagonal in 
Lorentz indices, i.e., is proportional to g^ u , which allowed us to suppress the Lorentz indices and write it in 
the form shown in Eq. (64). Such decomposition may not hold in other gauges. Inserting Eqs. (64) into the 
amplitude and defining 

M{q, A, k;p) = M\ + M 2 + M 3 + M 4 = (k) [Mf + M£ + M% + Mg] (66) 

where (k) is the polarization vector of the produced gluon. we get 

My = -gu(q) (t a S° F (q + k) r f (q + k,p) u(p) (67) 
M 2 = -gu(q)T f (q,p-k)S F (p-k)Jt a u{p) (68) 
M 3 = -gu(q) lv t h T b g a (k,p- q) u{p) G^(p - q) e„(fc) (69) 

M 4 = -g J dHu{q)T f { qi p-l)S%(p-l) lv t b T b g a (k,l)u(p)G^(l)e^k) (70) 

where r/ and T g are given by 

Tf (q,p) ee (27r)<5(p- - q~) 7 " / d 2 x t e *(«*-P*)-*« [V(x t ) - 1] (71) 

r g (q,p) ee 2p" (27r)5(p- - g") | d 2 ^ e'to"*)-*' [C/(x t ) - 1]. (72) 

The matrices V and f7 include all the multiple scatterings of the quark and gluon as they propagate in the 
strong classical field of the target and are given by 

V(x t ) = Pe i9 S dx ~ Ai(xt ''~ )ta (73) 
U(x t ) = Pe i9 I dz ~ Aiixt ' z ~ )T \ (74) 

Here t a and T a are matrices in the fundamental and adjoint representations of the SU (N) group respectively 
and A+(x t ,x-) = -gS(x-)^. 

With these definitions at hand, extracting an explicit factor of (27r)<5(p~ — q~ — k~) while using the delta 
function 5(l~ — k~) to do the l~ integration and performing the / + integration in M4 via contour integration 
using the (p — I) pole, we can write the amplitude as 
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M 2 
M 3 
M 4 



1 



'.</— -u{q)4{4 + #) 7" «(P) *° [^(«t + fct) - (2 7 r) 2 ( 5 2 (g t + fc t )] 

fc t )-(27r) 2 ^ 2 (9 t + fc t )] t a 
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"17^— T u (l) 7" ty-$)<{u{p) [V(q t 



2p-k 

W^- q u(<l) 1, u(p) d""(p - g) c M (fc) t b [*7 ba (<Zt + *t) - ^ (27r) 2 5 2 (g t + fc t )] 



(2tt 



-2 7" - /) 7„ --3— e M( fc ) 



P 

[K(«t + It) - (2^ 2 (g t + l t )]t b [U ba (h - l t ) - S ba {2nf5 2 (k t - h)] 

= _i 72 , 



(75) 



where d^ v (l) is related to the free gluon propagator via d^ v (l) = —i I 2 G^ u (I) and I — k , l + = — We 
have also set the transverse momentum of the incoming quark to zero without any loss of generality. We 
show the different diagrams contributing to the amplitude in Fig. 12. They correspond, respectively, to the 
quark multiply scattering from the target before or after radiating a gluon in Figs. 12-1 and 12-2 and the 
radiated gluon multiply scattering from the target in Fig. 12-3 while in Fig. 12-4 both the radiated gluon 
and the final state quark multiply scatter from the target. For the sake of clarity, momenta of the incoming 
quark and outgoing quark and gluon are shown explicitly in Fig. 12-1. 





(l) 



(2) 





(3) (4) 

FIG. 12. Diagrams corresponding to Mi, M 2 , M 3 , M 4 . 

In Fig. 13 we show one of the diagrams which are suppressed in the high energy limit and do not contribute 
and therefore, not included in Eqs. (75). This is due to the fact that the typical time scale for gluon emission 
is much longer than the time between rescatterings in the target nucleus. The diagram in Fig. 13 is, therefore, 
suppressed by a power of center of mass energy and can be safely neglected. Another diagram (not shown), 
which is suppressed in the high energy limit for the same reason, is when both the initial and final state 
quark lines as well as the radiated gluon multiply scatter from the target. 

To calculate the cross section, we need to square the amplitude \Ai\ 2 (66). There is a factor of — g^ + 
)[k/j,T)^ + ViJ.k u ] coming from squaring and summing over the polarization of the final state gluon in the 

light cone gauge which can be used to simplify the expressions. Furthermore, we define z = q~ /p~ so that 
1 — z = k~ jp~ . Below, we list the different contributions coming from squaring the amplitude. For reasons 
which will become clear shortly, we consider the square of Mi +M2 first. The contribution of the — g^ v term 
to the squared amplitude is (extracting a factor of g 1 for convenience) 



( J_ 

\T)-k 



-g^(M? + M^)t(Mf + M%) = I6p- P - 



z(l 



[zkt - (1 - z)q t f 
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+ 



z(l-z) 2 



+ 



Tr[Vt(g t + k t ) - (2n) 2 S 2 (q t + k t )]t a t a [V(q t + kt) - (2i:) 2 6 2 (q t + k t )\ 
Tr[V(q t + k t ) - (2t:) 2 5 2 (q t + k t )]t a t a [V\q t + k t ) - {2v) 2 5 2 (q t + k t )\ 

q 2 z 2 (l -z 2 ) l-z 



{1-zYd + z 2 )- 



+ 



k 2 \zk t -{\-z)q t f ' [zk t -(l-z)q t ] 2 
Trt a [V\q t + k t ) - (2ir) 2 6 2 (q t + k t )]t a [V(q t + k t ) - {2v) 2 5 2 {q t + k t )] 



(76) 



where Tr denotes trace of color matrices. This term is identical, up to color matrices, to the photon + 
quark production calculated in [35]. We now consider contribution of the [fc^r/j, + »7 M fe„] piece. Using the 
Dirac equation and the identity u(q)-y~u(p) — 2y / 2p~q~ (valid for on mass shell particles) simplifies the 
trace algebra considerably and we get 



M„ +VM (M ,, + M £)t( Ml " + Ml) = 32 p-p-z 2 

7] ' K 



V\qt + k t )t a -t a V\q t + k t ) 



w^^ tav{qt+kt) -h v{qt+kt)ta . 



(77) 



Note that this piece has no analog in QED and would vanish in the case of photon emission. 



^^^^^^^^ 

FIG. 13. A typical diagram which is suppressed in the high energy limit and, therefore, not included. 

A few remarks regarding the difference between photon and gluon radiation is in order here. In single 
inclusive photon production in pA collisions as considered in [35], the photon is emitted by a quark scattering, 
via multiple gluon exchanges, from a target which is treated as a classical gluon field, generated by recoilless 
sources of color charge. The photon current is conserved and satisfies /c M M M = due to gauge invariance. It 
is therefore enough to work in the covariant gauge where the sum over polarization of photons is just — g^ v . 
There is an essential difference between photon and gluon radiation here due to the fact that in the case of 
gluon radiation, one also needs to consider radiation of gluons from the target and not just from the quark. 
This is essential for gauge invariance of the amplitude and current conservation. However, it can be shown 
that as long as one works in the light cone gauge, the gluon radiation from the target vanishes identically. 
This is the case here since we are working in the light cone gauge. However, this means that one needs to 
keep the full projector — g^ v + (^)[fc jU r?i/ + Ty^fc^] rather than only the — piece which would be the case 
in the covariant gauge. We now consider the rest of the diagrams: 

\MtM,\ = Wp-p- z(l + Z 2 ) i~ Zqt ]\ qt + h ] 2 [U' ab (qt + k t ) S« b (2n) 2 6 2 (q t + k t )\ 

qt\zk t - (1 - z)q t y 

Trt b t a [V(q t + k t ) - (2ir) 2 6 2 (q t + k t )] 
\MlM 2 \ = \&p-p-z{\ + z 2 )^-[U^ b (q t + k t ) - S ab (27T) 2 S 2 (q t + k t )]Trt h [V{q t + k t ) - {2^) 2 5 2 {q t + k t )] t a 

|M 3 | 2 = 16p-p- Z{1+ 2 Z2) [U^ b (q t + k t ) - S ab (27T) 2 S 2 (q t + k t )}[U ca (qt + k t ) - 5 ca (2ir) 2 5 2 (q t + k t )] 
It 

Trt b t c 
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\M\M,\ = -I6p-p-z(l + z 2 ) J {^ ab (q t + h) - 5 ab (2n) 2 S 2 (q t + kt)] 



m\ma 



[U ca (k t - U) - 5 ca {2n) 2 5 2 {k t - l t )] Trt»{V(q t + l t ) - (2^ 2 5 2 {q t + l t )]t c 

-iep-p-z { i + z 2 ) J ^l^ttw ^ ~ h) - Sab{2 * mh ~ h)] 



Trt b [V\q t + It) - (27r) 2 5 2 (q t + l t )]t a [V{q t + k t ) - (2n) 2 5 2 (q t + k t )\ 
\M\M 2 \ = -16p-p-*(l + z 2 ) J ^y%^ [U ]aC {h - It) 6*%2n) 2 5 2 (k t - k)} 
Trt c [VHq t + It) - (2n) 2 S 2 (q t + k)][V(q t + h) - (2n) 2 5 2 (q t + h)]t a 
\M 4 \ 2 = 16p-p-*(l + z 2 ) J ^^ l j^[U^(k t - It) S^(2n) 2 6 2 (k t - l t )} 
[U ab (k t - l t ) - S ab (27r) 2 S 2 (k t - l t )]Trt c t b [V\q t + l t ) - (27r) 2 S 2 (q t + l t )] 

[V(qt + lt)-(2ir) 2 5 2 (q t + lt)]. (78) 

Note that for a given interference term such as l-MgMiJ, there is also the conjugate term |M^M 3 | which 
is obtained from |M^Mi| by daggering the color matrices. Eqs. (76-78) provide the complete expression 
for the amplitude squared \M\ 2 . In order to get the invariant cross section, one needs to include the phase 
space and the flux factors given by ^tt) 3 2cr ' (2C) 3 2F 7, ' an< ^ 2p^- Including a factor of 1/2 coming from 
averaging over the incoming quark spin and restoring the coupling constant and the overall delta function, 
the invariant cross section is given by 

da qA ^ qgX 1 1 

« k ~^q-¥k- ^ I<F W {2W)S(P ~ ~ " ~ k ~ )92 |M|2 ' (79) 

This is the invariant cross section for production of a quark and gluon in the scattering of a quark on a 
target nucleus (or a proton at small x) including classical multiple scattering. In order to get the invariant 
cross section for production of two hadrons or two jets in a proton (deuteron)-nucleus collisions, one needs 
to convolute the cross section given in Eq. (79) with the (valence) quark distribution function of a proton 
or deuteron and the quark or gluon fragmentation functions 

da pA ^ hl h2 x da qA ^ qgX 
EhA d^k h2 = *W 8 «" fc "^-®^W ^ (80) 

where q P {x q ) is the quark distribution function in a proton and the quark and gluon fragmentation functions 
are denoted by D^^zi) and D g l2 (z2), while ® denotes a convolution over Bjorkcn x for distribution function 
and over z\,z 2 for fragmentation functions. The cross section calculated here is valid when one produces 
two hadrons (or jets) in the forward rapidity region of a proton (deuteron) nucleus collision. It includes the 
effects of quantum evolution (in x) in the target. To see this one has to evaluate the Wilson lines (J7's and 
V's) correlators in Eqs. (76-78) using the JIMWLK evolution equation [10]. One can use this cross section in 
order to investigate two particle correlations (back-to-back jets) in the RHIC forward rapidity region which 
can be measured, at RHIC for example, by the STAR detector at rapidity y = 3.8. 

In order to consider the case when one of the hadrons in produced in the mid rapidity region, one needs 
to allow the possibility that the gluon is radiated not from the valence quark directly, but from anywhere 
along the (in principle, non-linear ) gluon cascade between the valence quark and the target. 

At this point, it is worthwhile to make a connection between the notations used in different sections since 
they may seem disjoint to a casual reader. The degrees of freedom are indeed the same even though they 
are denoted differently due to convenience. In Sections II and III the forward scattering amplitude of a 
quark-anti-quark dipole on the target is denoted by N(x t ,yt) and can be expressed as 

N(x u y t ) = 1 - ^ Tr(vHx t )V(yt)) , (81) 

where V is a path-ordered integral in the fundamental representation, used in section TV and defined in 
Eq. (73). Also, the adjoint dipole amplitude denoted Nq in, for instance, Eq. (21), is equal to 
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N G (x t ,y t ) = l- J J-^Tr(uHx t )U(y t )), 



(82) 



where U is the path-ordered integral in the adjoint representation, used in section IV and defined in 
Eq. (74). Furthermore, the S'-matrix of the color quadrupole interaction with the target, which is de- 
noted Qo(xt,Ut, z t ,r t ) in the classical case and calculated in Eq. (14) and denoted Q(x t ,yt, z t ,r t ) in the case 
of quantum evolution included in Eq. (28), can be rewritten in terms of the correlator of four path-ordered 
exponentials 

Q(x t ,y t ,z t ,r t ) = -L Tr (v\x t )V{z t )V\r t )V{y t )) . (83) 

The relations in Eqs. (81), (82) and (83) between N, N G , Q and the correlators of V's and U's hold even 
when the quantum evolution is included. 



V. CONCLUSIONS 



In this paper we have calculated two cross sections for inclusive two-particle production relevant for the 
dAu run at RHIC and for the upcoming pA run at LHC. The cross section for two-gluon production at mid- 
rapidity for DIS is given by Eq. (32). The expression in Eq. (32) includes all multiple rescatterings of the 
produced gluons on the target, along with the non-linear small- a; evolution effects [9]. Even though, unlike 
the single gluon inclusive production cross section of [23] [see Eq. (19) above], our two-gluon cross section 
in Eq. (32) can not be cast in fcT-factorized form, it can be easily generalized to the case of proton-nucleus 
(pA) scattering. Following [23] we note that the probability ni(x Q, x 12 , b, y) of finding a dipole 12 in the 
original dipole 00 can be related to the unintegrated gluon distribution <f>(q, y) as 

f dH-^n^x^y) = ^ / rfV^^fp), (84) 

where the coefficient on the right hand side of Eq. (84) has been fixed in order for Eq. (19) to be reducible to 
conventional fey-factorization form of [37] . The information about the original dipole 00 is now contained in 
its unintegrated gluon distribution 4>{q, y). Eq. (84) makes generalization of Eq. (19) from DIS to pA rather 
straightforward: instead of the unintegrated gluon distribution function <f>(q, y) of the incoming dipole, 
one has to use a BFKL evolved unintegrated gluon distribution <p(q,y) of the proton in Eq. (84), and, 
consequently, in Eq. (19). 

To repeat the above procedure for Eq. (32) we have to devise a generalization procedure for the probability 
of finding two dipoles n 2 as well. To do that, let us first clarify the physical meaning of Eq. (84). The 
term in Eq. (84) is due to V^ 12 , which is usual to the definition of unintegrated gluon distribution in terms of 
the dipole amplitude (see Eq. (2) in [2]), and 1/V* 12 , which is proportional to gluons' propagators in a two- 
gluon exchange amplitude. Thus, in Eq. (84) the gluon distribution is obtained from the dipole probability 
ni by connecting two i-channel exchange gluons to the dipole 12 in it. Now, generalization of Eq. (84) to n 2 
becomes manifest: one has to connect two exchange gluons to each of the two produced dipoles. The final 
expression reads 



l d b 22 ^ 2 ^ 2 n2(x J0 , x_ Q, Y; x l7 x_ j, y\, x 2 , x_ 2l 1/2) 



(2tt 



'11 • L 22 



J d 2 qd 2 le i ^ +a --^4>2(g,Y-yv,LY-y 2 ), (85) 

where 6 x j = {x 1 +x_i)/2, b 22 = (x 2 +x 2 )/2, and (f> 2 (q,Y — y\;l,Y — j/2) is the two-gluon distribution function 

in the incoming dipole 00, with the two gluons having transverse momenta q and I and rapidities Y — y\ and 
Y — y 2 with respect to the projectile onium. 

Analyzing Eq. (32) one can see that for scattering on a large nucleus both m and n 2 come into Eq. (32) 
integrated over impact parameter(s), as employed in Eqs. (84) and (85). Therefore, using Eqs. (84) and (85) 



27 



one can rewrite Eq. (32) in terms of single and double unintegrated gluon distributions <j> and <p 2 . Taking 
these distributions for a proton (deuteron) instead of the quarkonium would accomplish generalization of 
Eq. (32) to the case of p(d)A scattering. 

Eqs. (76-78) along with Eq. (79) give us a production cross section for a valence quark and a gluon in the 
forward rapidity direction in p(d)A scattering. If the correlators of Wilson lines in Eqs. (76-78) are averaged 
in the Gaussian approximation [8], the obtained cross section (79) would reduce to the quasi-classical result 
containing multiple rescatterings only. (For an explicit evaluation of color averaging of the Wilson lines using 
a Gaussian weight, see [21,57].) If the Wilson lines are averaged with the weight function obtained from 
solving the JIMWLK evolution equation [10], than Eq. (79) would include the complete effects of small- a; 
evolution as well. 

Before we conclude we would like to make a comment about the applicability of Eqs. (18), (76-78) for RHIC 
kinematics. Indeed, in deriving these equations, we have assumed for simplicity that the gluons are widely 
separated in rapidity, y 2 ^ j/i • On the other hand we know that particle production at mid-rapidity at RHIC 
appears to be better described by the quasi-classical physics leading to Cronin enhancement. Therefore, if 
one of the produced particles is at forward rapidity with the other one being at mid-rapidity our formulas 
would apply, though one would not need to include the small-x evolution between the target nucleus and 
the particle produced at mid-rapidity, since there the physics is quasi-classical. However, the suppression in 
R dAu : which is most likely caused by small-x evolution, sets in already at rapidity r/ = 1 and continues all 
the way up to the highest achievable rapidity at RHIC [12,13,15]. That means quantum evolution describes 
physics at r\ > 1. Therefore, if both of the produced particles are at rapidity r\ > 1, say if y\ = 1 and y 2 = 3, 
we can still have a large rapidity interval between them, y 2 ^> yi, and have quantum evolution between the 
target and the gluon at yi included in Eqs. (18), (76-78). Indeed, the upcoming pA run at the LHC would 
have a much wider rapidity window, where our results would be even more applicable. 
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APPENDIX A 

In this Appendix we calculate the S'-matrix of the interaction a quadrupolc 2,2', 1,1' with the target 
nucleus. The S'-matrix includes the Glauber-Mueller [44,7,17,48] multiple rescatterings only and is denoted 
by Qo {x 2 ,x 2 i, m Sect. II above. The possible interactions are shown in Fig. 14. The first term there 
corresponds to the case where all the interactions in the amplitude and in the complex conjugate amplitude 
are virtual, i.e., each nucleon exchanges two gluons with the qq pair and remains intact. This is the diffractive 
piece of the interaction [48]. The gluons connect to both the quark and the antiquark lines. This is denoted 
by leaving the gluons lines disconnected at the top ends. All the virtual exchanges are leading at large N c . 
The first diagram in Fig. 14 gives a contribution 

e -x 2 21 ln(l/x 21 A) Q 2 s0 /A e -x 2 2ll , \n(l/x 2/1 , A) Q 2 s0 /4 

which is just a product of the S-matrices of the dipoles 12 and 1'2'. 

The second diagram in Fig. 14 corresponds to the case of at least one real interaction: there the nucleon 
at longitudinal coordinate z interacts with the qq pair by a single gluon exchange in the amplitude and in 
the complex conjugate amplitude. The single gluon exchange breaks up the nucleon in the final state. We 
will refer to this interaction as real [17,20]. The interaction of the nucleon at z is chosen to be the first real 
interaction: all prior exchanges are virtual (exchanges to the left of z in the amplitude and to the right of z 
in the complex conjugate amplitude). After the interaction of nucleon at z the exchanges can be both real 
and virtual. However, in the large N c limit, only those real exchanges contribute where gluons connect to 
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either lines 1 and 1' or lines 2 and 2'. The color structure is similar to the dipole model [49]: the nucleon at 
z splits the original single quark loop 2, 2', 1, 1' into two, and the successive interactions can only take place 
within each of the two resulting loops. 




FIG. 14. Leading diagrams contributing to the interaction of a color quadmpole 2,2', 1, l' with a nuclear target in 
the large- N c limit. 



The contribution of the second graph in Fig. 14 is therefore 

^ *1 e -i MV^iAHx^, in(i/* a 'i'A)] Q 2 s0 i * [ x l 2 , ln(l/x 2 2'A) + x\ v In(l/a: U /A) 



+a| x , ln(l/x2i'A) + a&i ln(l/a; 2 /iA)] e~* ln ( 1 /^ 2 'A)+^ 1 , in(V*u'A)] Q= V. (A2) 

In Eq. (A2) the first exponent resums all virtual interactions before the first real interaction at z, the 
second exponent resums all the real and virtual interactions with dipoles 22' and 11' following the first real 
interaction, and the term in between accounts for the first real interaction itself. We also average over the 
longitudinal coordinate z which varies from to L, where L is the longitudinal extent of the nucleus at a 
given impact parameter. 

Performing the integration over z in Eq. (A2) and adding to it the contribution from Eq. (Al) yields 
Eq. (14) in the text. 
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